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ABSTRACT 

We present a quantitative and homogeneous analysis of the broadband multi- 
color photometric data sets gathered so far on rapidly pulsating hot B subdwarf 
stars. This concerns seven distinct data sets related to six different stars. Our 
analysis is carried out within the theoretical framework developed by Randall et 
al., which includes full nonadiabatic effects. The goal of this analysis is partial 
mode identification, i.e., the determination of the degree index / of each of the 
observed pulsation modes. We assume possible values of / from to 5 in our 
calculations. For each target star, we compute a specific model atmosphere and 
a specific pulsation model using estimates of the atmospheric parameters coming 
from time-averaged optical spectroscopy. For every assumed value of /, we use a 
formal approach to model the observed amplitude-wavelength distribution of 
each mode, and we compute a quality-of-fit Q probabihty to quantify the derived 
fit and to discriminate objectively between the various solutions. We find that no 
completely convincing and unambiguous / identification is possible on the basis 
of the available data, although partial mode discrimination has been reached for 
25 out of the 41 modes studied. A brief statistical study of these results suggests 
that a majority of the modes must have / values of 0, 1, and 2, but also that 
modes with / = 4 could very well be present while modes with / = 3 appear to 
be rarer. This is in line with recent results showing that / = 4 modes in rapidly 
pulsating B subdwarfs have a higher visibility in the optical domain than modes 
with / = 3. Although somewhat disappointing in terms of mode discrimination. 
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our results still suggest that the full potential of multicolor photometry for I iden- 
tification in pulsating subdwarfs is within reach. It will be a matter of gathering 
higher S/N ratio observations than has been done up to now. 

Subject headings: stars: stars: horizontal-branch — stars: interiors — stars: 
oscillations — subdwarfs 



1. INTRODUCTION 

It is well established that an analysis of multicolor photometric data, whereby one 
compares oscillation amplitudes and/or phase differences in two or more wavebands, can 
lead to partial mode identification in pulsating stars (see, e.g., Heynderickx, Waelkens, & 
Smeyers 1994). The technique has been used successfully in the past to infer the values 
of the degree index / of pulsation modes for many types of oscillating stars. To name just 
a few, let us mention that it has been applied to 5 Scuti stars (Garrido, Garcia-Lobo, & 
Rodriguez 1990), P Cepheids (Cugier, Dziembowski, & Pamyatnykh 1994), ZZ Ceti white 
dwarfs (Fontaine et al. 1996), 7 Doradus stars (Breger et al. 1997), and Slowly Pulsating 
main sequence B stars (Dupret et al. 2003). Attempts to understand the observed period 
spectrum on the basis of multicolor photometry have also been made for EC 14026 stars, 
first by Kocn (1998) and more recently by Jcffcry ct al. (2004; 2005). The former study gives 
a qualitative interpretation of the periods observed for KPD 2109+4401 based on the theory 
of Watson (1988), while the latter assert to have provisionally identified or constrained 
the I values for some of the modes detected for the fast oscillators KPD 2109-1-4401, HS 
0039-1-4302, and PG 0014-1-067. In these cases, the I identification is based on a qualitative 
comparison with amplitude ratios computed by Ramachandran, Jeffery, & Townsend (2004) 
in the adiabatic approximation for representative models. 

A theoretical framework for the quantitative exploitation of multicolor photometry for 
pulsating sdB stars has been put forward recently by some of us (Randall et al. 2005). In 
that paper, the potential of the technique as applied to both types of pulsating sdB stars 
(the short-period p-modc pulsators of the EC 14026 type and the slowly oscillating ^f-mode 
variables of the PG 1716 type) has been explored in detail. The method of Randall et al. 
(2005) features a full nonadiabatic treatment of the atmospheric layers and uses a designated 
model atmosphere code which automatically incorporates the wavelength dependence of the 
limb darkening — thus avoiding the need for approximate parameterized limb darkening 
coefficients as used in most other multicolor photometric studies with the notable exception 
of Ramachandran et al. (2004). Of central importance, the method of Randall et al. (2005) 
has led to the first (and, so far, only) unambiguous determination of the I index of a pulsation 
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mode in a oscillating sdB star using multicolor photometry, thus demonstrating the feasibility 
of the approach. 

Indeed, at the time of the writing of the Randall et al. paper, the best available am- 
phtude estimates of oscillation modes in a pulsating EC 14026 star were those based on the 
superb u'g'r' data set gathered by Jeffery et al. (2004) on KPD 2109+4401 {V = 13.38) 
with t7LTJ?A CAM mounted on the 4.2-m William Herschel Telescope (WHT). Considering 
the largest amplitude mode reported with a period of 182.42 s, and adopting at face value 
the three amplitudes and their quoted uncertainties {u' = 8.87 ±0.04 mmag, g' = 6.58 ±0.04 
mmag, r' = 6.15 ± 0.04 mmag), Randall et al. (2005) were able to demonstrate that this 
mode must be a radial mode (/ = 0) as first indicated by Jeffery et al. (2004). 

The quality of these unique observations has remained unsurpassed so far, but there 
are nevertheless other available multicolor data sets that certainly deserve to be analyzed 
in the same fashion now that there is a proven theoretical framework to do that. While 
it is known that it is generally difficult to discriminate between modes with I — 0, 1, or 
2 on the basis of multicolor photometry for EC 14026 stars (Ramachandran et al. 2004; 
Randall et al. 2005), one may hope to distinguish between modes of that group from modes 
with Z = 3 or / = 4 or possibly / = 5. Our survey of the published work in the field has 
revealed the following available data sets 1) the UBVR photometry of Koen (1998; herafter 
K98) on KPD 2109±4401, 2) the UBV photometry of Silvotti et al. (2000; herafter SOO) 
on PG 1628+563B, 3) the BUSCA photometry of Falter et al. (2003; hereafter F03) on 
PG 1605±072, 4) the u'g'r' / ULTRAC AM photometry of Jeffery et al. (2004; hereafter 
J04) on HS 0039±4309, 5) the u'g'r' / ULTRAC AM photometij of Jeffery et al. (2004) on 
KPD 2109+4401, 6) the Stromgren photometry of Oreiro et al. (2005; hereafter O05) on 
BAL 090100001, 7) the UBVR photometry of Baran et al. (2005; hereafter B05) on BAL 
090100001, 8) the -u'^V'/ ULTHACAM photometry of Jeffery et al. (2005; hereafter J05) on 
PG 0014+067, and 9) the w'^V/ f/LTRACAM photometry of Aerts et al. (2006) on SDSS 
J171722. 08+58055. 8. Except for cases 5) and 9) (see below), we analyze all of the other data 
sets in the present paper. 

Partial mode identification (the identification of the degree index I) through multicolor 
photometry is a worthwhile venture in itself, but it has become particularly important for 
EC 14026 pulsators now that specific asteroseismic models have been proposed for a few of 
them. Indeed, by combining the forward method with high S/N ratio spectroscopy, it has 
been possible to carry out complete asteroseismological analyses for four EC 14026 pulsators 
so far: PG 0014+067 (Brassard et al. 2001), PG 1047+003 (Charpinet et al. 2003), PG 
1219+534 (Charpinet et al. 2005a), and Feige 48 (Charpinet et al. 2005b). The ultimate 
product of these analyses is the determination of the global structural properties of the stars. 
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but complete mode identification (i.e., a determination of the radial order k and the degree 
index / of each pulsation mode) is a byproduct of the method. Multicolor photometry 
can thus be used as testing grounds for the proposed seismic models by checking if the / 
assignments are correct. Of the four EC 14026 pulsators with a proposed specific seismic 
model, only PG 0014+067 has been observed so far in several filters simultaneoulsy (J05). 
We will discuss this interesting case below as part of our analysis of that data set. 

The main purpose of the paper is to present a quantitative and homogeneous analysis 
based on the method of Randall ct al. (2005) for each of the available multicolor data sets that 
we found in the published literature on pulsating EC 14026 stars. Considerable observational 
efforts went into the gathering of these data sets and we feel that it is most worthwhile to 
attempt extracting the maximum information from them. Partial mode identification is 
the goal. We note that the exact same approach has been used recently by Aerts et al. 
(2006) for analyzing their ULTRACAM/WRT data of the very faint EC 14026 pulsator 
SDSS J171722. 08+58055. 8 {B ~ 16.7), so it is not necessary to reconsider their data here. 
Unfortunately, in that case, the identification of I was not possible for the two pulsation 
modes uncovered because the uncertainties on the observed amplitudes turned out to be 
too large due to insufficient S/N ratio. In contrast, some interesting mode identifications 
appear possible for the Jeffery et al. (2004) data set on KPD 2109+4401 (beyond that for 
the 182.42 s mode discussed in Randall et al. 2005), but we defer our analysis to a separate 
publication that will include unpublished additional observations that we gathered at the 
Canada-Prance-Hawaii 3.6-m Telescope on that star. We note also that two-color (U and 
R) photometry has been gathered and analyzed for two PG 1716 sdB pulsators in order to 
constrain the / index (sec the paper of Randall et al. 2006a on PG 1627+017 and that of 
Randall et al. 2006b on PG 1338+481). 



2. METHOD 

2.1. The computations of theoretical ampHtudes 

To analyze multicolor data for pulsating sdB stars, we follow the modeling method 
developed by Randall et al. (2005). While we refer the reader to that paper for more 
information, we include here a brief recap of the most relevant aspects for convenience. The 
technique employed is based on the well-known fact that the observed wavelength dependence 
of an oscillation's amplitude and phase bears the signature of the mode's degree index I 
and its period, as well as the star's atmospheric parameters, the intrinsic amplitude of the 
periodicity, and its viewing aspect. To first order, the influence of the latter two (unknown) 
parameters can be eliminated by computing the ratio of pulsational amplitudes (and the 
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difference between phases) as measured in two or more bandpasses. Given a target's effective 
temperature and surface gravity, the observed mode's degree index can thus, in principle, be 
inferred from multicolor photometry, leading to partial mode identification. 

The method proposed by Randall et al. (2005) incorporates a full nonadiabatic descrip- 
tion of the atmospheric layers in the computations of theoretical pulsation observables. As 
in main sequence stars (e.g., Dupret et al. 2003), such description is found to be quite im- 
portant in sdB stars as well. In particular, while the predicted phase shifts between various 
bandpasses generally remain quite small — they would all be identical to zero in the adia- 
batic approximation — Randall et al. (2005) found that amplitude ratios cannot, in that 
approximation, be computed with enough accuracy for quantitative studies. Furthermore, 
they found that the amplitude ratios do depend sensitively on the atmospheric parameters of 
the target, so for the purposes of our present analysis, we need to compute a specific detailed 
model atmosphere as well as a full pulsation model for each one of our targets. 

For each value of I, the brightness variation expected across the visible disk during 
a pulsation cycle can be expressed in terms of temperature, radius, and surface gravity 
perturbations to the emergent flux. These in turn are dependent on quantities obtainable 
from model atmospheres and nonadiabatic pulsation theory, as well as on the period of the 
mode in question (see, e.g., equations 34 — 38 of Randall et al. 2005). The model atmosphere 
parameters are made up of the logarithmic derivative of the emergent flux with respect to 
the effective temperature / surface gravity and the weighted limb darkening integral together 
with its derivatives. Initially computed as monochromatic quantities {aru, (^gu, i>iui biu,T-, and 
hiy^g in the notation of Randall et al. 2005) from a specially modified sdB atmosphere code, 
they are subsequently integrated over the bandpasses of interest (the quantities axx) ctgx, 
bix, bix^T, and bix,g, where x symbolizes the appropriate bandpass, again in the notation of 
Randall et al. 2005) to allow for comparison with observations. 

The second set of parameters includes the nonadiabatic quantities R and ^^t, which are 
related to an oscillation's departure from adiabacity in amplitude and phase respectively. 
Together with the adiabatic gradient, Vad, they are computed on the basis of the so-called 
"second generation" static stellar models described in Charpinet et al. (2001). These full 
models are submitted to adiabatic and nonadiabatic pulsation calculations that yield, among 
other things, the relative behaviour of the radius (gravity) and temperature perturbations 
in the stellar atmosphere for each mode considered. Integrated over the atmospheric layers 
contributing most to the emergent flux (taken to lie at optical depths r = 0.1 — 10), the ratio 
of the two perturbations' moduli and their phase lag yield the atmosphere-averaged values 
of R and describing the departure from adiabacity of the observed brightness variations 
for each mode. The atmospheric value of Vad is obtained in a similar way. 
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For each reported mode in each observed target, we calculated the predicted amplitude 
ratios and phase differences from the available bandpasses for values of the degree index from 
/ = to / = 5. Since, as it turns out, there is very little information content in the phase 
shifts (the predicted differences are indeed too small in all cases of interest to be detectable 
at the accuracy achieved in the available data sets), wc concentrate solely in what follows 
on comparing predicted and observed amplitudes. To provide a quantitative framework for 
such a comparison, we contrast for each mode the predicted multicolor amplitudes with 
those observed using a minimization procedure following Fontaine et al. (1996). For 
every degree index I, the theoretical amplitudes alf^^g in each of the A'" available bandpasses 
i are multiplied by a scale factor fi chosen in such a way as to minimize 

1=1 

where a^^^ is the amplitude observed in a given waveband and cr* is the error on the measure- 
ment. Compared to the standard normalization of all amplitudes to one particular waveband, 
this is a more objective way of determining the quality of a match, since the data from all 
bandpasses are treated on the same footing. 

Finally, since our procedure is a standard approach, it is possible, following Fontaine 
et al. (1996), to quantify the value of a given solution by computing explicitly the 
quality-of-fit Q as described in Press et al. (1986). The quantity Q depends on the value of 
for each solution and the number of degrees of freedom (A^ fitted points minus the free 
parameter j] gives — 1 degrees of freedom in the present case). We adopt the canonical 
notion suggested by Press et al. (1986) that a fit is acceptable if its quality-of-fit Q > 0.001. 
For a system with a single degree of freedom, this value of Q corresponds to the case where 
the predicted amplitudes fall simultaneoulsy in all bandpasses of interest within 3(t\ Hence, 
we shall use the criterion Q > 0.001 to determine if a given x^ fit is acceptable and to 
discriminate quantitatively betweeen possible I solutions. 



2.2. The available data sets 

The vital characteristics of each of the seven data sets used in this study are summarized 
in Table 1. The first column gives the name of the EC 14026 pulsator, the second one gives 
its V magnitude {B if V is unavailable), the third one hsts the total length of the data set, 
the fourth column refers to the bandpasses of interest, the fifth one indicates the instrument 
used, and this is followed by the identification of the telescope used and its location, and 
finally a reference to the original paper is given. 
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In the calculations of theoretical amplitudes for each target star, it is necessary to specify 
the atmospheric parameters of the star, log g and Tgfj. To do this, wc rely on the recent 
and ongoing spectroscopic efforts of Green, Fontaine, & Chayer (2006, in preparation) to 
derive reliable and homogeneous atmospheric parameters for a large sample of sdB stars. 
For each of our target stars, we used estimates for the atmospheric parameters coming from 
that study. These are given in Table 2. 

In addition, we indicate in the same table the three ingredients that we used in the con- 
volution process transforming monochromatic quantities into waveband- integrated quantities 
suitable for comparison with the observational data. These include the transmission curve of 
each filter, the response curve of the detector, and the extinction curve of the site. Tests in- 
dicate that, by far, the transmission curves of the filters are the most important component, 
and that, at the other extreme, the choice of the extinction curve is a sophistication that 
has little impact for sites located at roughly the same altitude. In this connection, we picked 
the atmospheric transparency curve of a representative medium altitude site (in this case, 
Kitt Peak National Observatory), a choice that should be appropriate for all but perhaps 
Mt. Suhora Observatory in Poland located at an altitude of about 1000 m. Further tests 
indicate that the largest sources of uncertainty for the computed amplitudes are associated 
with the uncertainties on T^s and log g, assuming that the filters used by the observers have 
indeed the claimed transmission properties. For typical uncertainties on the atmospheric pa- 
rameters of ATeff = ± 400 K and Alogg' = ± 0.05 (this excludes possible systematic effects 
associated with the model atmospheres used by Green et al. 2006), the real effects on the 
theoretical amplitudes remain quite small and generally not significant enough to affect the 
implications for mode identification (given the observational uncertainties on the observed 
amplitudes) . 

For each data set, we generally accepted the data, i.e., the observed multicolor am- 
plitudes and their associated uncertainties at face value. This implicitly implies that each 
reported mode was fully resolved and that its amplitudes were free of spectral contamination 
from neighborhing peaks in the Fourier domain. This is an important issue: if a mode is not 
sufficiently "monochromatic" in frequency, attempts to model its / signature as a function 
of wavelength are usually frustrated. One exception to this is the case of an unresolved 
frequency multiplet due to rotational splitting. In that case, since the amplitudes do not de- 
pend on the azimuthal order m, the modeling remains valid for such an unresolved multiplet 
in a slowly rotating star. 

An equally important concern is the question of the uncertainties on the observed am- 
phtudes. Since our approach is based on a statistics, it is absolutely essential to have 
realistic and accurate estimates of these uncertainties. The quality of the fit, as measured by 
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the quantity Q, and which gives the discriminatory power of the method, depends directly 
on the values of these uncertainties. We found, in this connection, that our failure sometimes 
to model in an acceptable way a given pulsation mode could be traced back to the fact that 
the errors on the amplitudes were probably underestimated. 

3. RESULTS 

3.1. KPD 2109+4401 (K98) 

The paper of K98 is the first one reporting on multicolor photometry for a rapidly 
pulsating sdB star. It was a key source of motivation and inspiration for the theoretical work 
that we presented in Randall et al. (2005). In K98, the author reported, among other things, 
on some 27 h of UBVR photometry on KPD 2109+4401 {V = 13.38) gathered with a 4- 
channel photometer mounted on the McDonald Observatory's 2.1-m Struve Telescope. Seven 
distinct pulsation modes were uncovered, and U BVR amplitudes and phases were provided 
for each one of them. The author, however, was clearly reluctant to provide estimates 
of the uncertainties on the amplitudes and phases, being quite aware, as he argued, that 
formal errors coming from least-squares fits of light curves tend to be unreliable. Reading 
between the lines, and using the information on the relative efficiency of each bandpass of the 
Steining photometer that we had from previous work (Fontaine et al. 1996), we established 
that reasonable estimates of the observational uncertainties on the amphtudes would be 
cr{U) — 0.105 mmag, cr(-B) = 0.105 mmag, (t{V) — 0.085 mmag, and cr{R) — 0.155 mmag, 
irrespective of the amplitude of a given mode. 

Table 3 summarizes the results of our modeling effort for the K98 data sct.^ For each 
observed mode, the table gives a block of data. The first fine in each block gives the period 
of the mode of interest, and this is followed by the observed amplitude and its uncertainty 
in each of the available bandpass. The next fine gives the best-fit theoretical amplitude in 
each of the bandpass for an assumed degree index / = 0, and this is followed by the value 
of obtained in the minimization procedure (eq. 1), and the value of the quality-of-fit Q. 
The next five lines show similar data, but for assumed values of / from 1 to 5, respectively. 

Figure 1 is a graphical representation of our results for the 182.42 s mode (the first 
one fisted in Table 3 and the second largest ampfitude one in the K98 data) . The observed 



^Note that we explicitly used the transmission curves of the UBVR filters of the Steining photometer 
as published in Robinson et al. (1995) for this analysis. These are similar but not exactly the same as the 
standard Johnson-Cousins filters. 
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behavior of that mode is well modeled for values of Z = 0, Z = 1, or / = 2 as can be seen in the 
figure. The plot is consistent, of course, with the values of Q all larger than 0.001 for those 
three possibilities, while it also clearly shows the poor fits {Q <^ 0.001; see Table 3) obtained 
for models assuming / = 3, / = 4, or / = 5. We thus find that, while the uncertainties on 
the observed amplitudes are too large for us to discriminate between the values of / from 
to 2, we can safely discard the possibilities that the 182.42 s mode in KPD 2109+4401 has 
a degree index of 3, 4, or 5. We note that this result is consistent with the determination of 
/ = proposed by Randall et al. (2005) for that mode on the basis of the the much higher 
S/N ratio data of Jeffery et al. (2004). 

The highest amplitude mode in the K98 data is the 196.31 s pulsation which, interest- 
ingly enough but not uncommonly in sdB pulsators, showed smaller amplitudes than the 
182.42 s mode when KPD 2109+4401 was observed by Jeffery et al. (2004). Our modeling 
of the K98 data shows a very similar situation to the previous mode as can be seen in Figure 
2. Again, we cannot formally discriminate between / = 0, / = 1, or / = 2, but we can safely 
exclude the possibilities that the 196.31 s mode in KPD 2109+4401 has a value of / = 3, 
/ = 4, or / = 5. The very small values of Q associated with these solutions (Table 3) certainly 
confirm this conclusion. 

The most interesting case in the K98 data is that of the 198.19 s mode (the third 
largest amplitude mode), which shows an amplitude- wavelength behavior characteristic of 
a / = 4 mode as can be appreciated in Figure 3. In fact, formally speaking, the / = 4 
solution is the only one acceptable, with a value of Q > 0.001 in Table 3. While this 
remains highly suggestive, we caution that one should not jump too hastily to the conclusion 
that mode discrimination has been achieved beyond any doubt for that mode. The reason 
for our cautionary remark is that the next best fit, the one corresponding to / = 2 and 
Q — 2.77 X 10~^ in Table 3, does not show a quahty-of-fit value that is much smaller than the 
passage criterion of Q — 0.001. For instance, if we were to increase the uncertainties on the 
observed amplitudes by some 10%, the solution with / = 2 would become formally acceptable. 
So here is a case where realistic and accurate estimates of the amplitude uncertainties become 
critical. 

For the other modes uncovered in K98, we formally constrain the value of I to 0, 1, 2, 
or 4 for the 184.72 s pulsation, to 2 or 4 for the 184.75 s pulsation, and to 1, 2, or 4 for the 
191.85 s pulsation. No mode discrimination was possible for the lowest amplitude mode at 
196.69 s, and all values of / from to 5 provide acceptable fits. We note that the provisionary 
/ assignments made by Jeffery et al. (2004) on the basis of their ULTRACAM/WYIT obser- 
vations are consistent with our more quantitative results here for all modes, except perhaps 
for the 198.19 s pulsation for which we may have isolated a unique value at Z = 4, although. 
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as discussed just above, the I — 2 solution (which is the suggestion of Jeffery et al. 2004 for 
that mode) should not be dismissed too hastily. 

3.2. PG 1618+563B (SOO) 

Amplitude and phase data were reported for two modes (139.3 s and 143.9 s) observed 
in PG 1618+563B {V = 13.52) by SOO. Those were derived from UBV observations using 
the 3-channel Tromso photometer attached to the 2.5-m NOT Telescope. Unfortunately, 
only 1.5 h of data were gathered and, consequently, the S/N ratio was not large enough to 
allow / index discrimination. Our results are summarized in Table 4 (the format is similar 
to that of Table 3) and clearly show that all possible / values from to 5 provide acceptable 
fits for both modes. Figure 4 gives a graphical representation of our results for the 139.3 s 
mode. 



3.3. PG 1605+072 (F03) 

The relatively bright {V = 12.92), large-amplitude sdB pulsator PG 1605+072 has been 
observed by F03 who provided 12.3 h of BUSCA photometry, a 4-channel system allowing the 
simultaneous observations in bandpasses centered on 3600 A{UV), 4800 A{B), 6300 A(i?), 
and 8000 A(NIR). These data were gathered at the 2.2-m telescope at Calar Alto. F03 
isolated 11 modes in their light curves and provided measurements of amplitudes and phases 
in their Tables 2 and 3. However, only 5 of those modes have simultaneous measurements 
in all the 4 BUSCA bandpasses. We have restricted our analysis to these 5 modes. 

In our modeling effort, we were careful to convolve the transmission curves of the 4 
filters with the response curve of the BUSCA instrument as provided to us by O. Cordes 
(2005, private communication). We also point out that the estimates of the atmospheric 
parameters that we used for PG 1605+072 (see Table 2) are remarkably consistent with the 
independent values obtained by Heber, Reid, & Werner (1999), which are Tes — 32,300 ± 
300 K and log g = 5.25 ± 0.05 dex. 

Despite careful modeling, we were unable to reproduce satisfactorily the amplitude- 
wavelength behaviors of all 5 modes when using the amplitudes and their uncertainties given 
in Tables 2 and 3 of F03. In all 5 cases, and for all assumed values of I from to 5, 
we obtained unacceptable quality-of-fit values of Q <^ 0.001. While we cannot completely 
rule out the possibility that our modeling is inadequate for this particular star, we believe 
instead that this problem is related to the fact that the amplitude uncertainties quoted 
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by F03 most likely underestimate the true errors. In fact, the authors of F03 seem to be 
aware of that and provide two other different estimates of the amplitude uncertainties. For 
instance, they explicitly write that the quoted uncertainties on the reported amplitudes in 
their Tables 2 and 3 are formal fit errors from a least-squares sine fit procedure, while the 
(much larger) uncertainties shown in their Figure 5 for 4 different modes have been estimated 
in another way. Furthermore, F03 write in their subsection 3.2 that the la noise level in 
their photometry corresponds to 1.52 mmag in the UV bandpass, 1.53 mmag in B, 1.12 
mmag in R, and 1.37 mmag in NIR, which is clearly at odds with the much smaller formal 
fit errors reported in their tables. 

We used these last figures to derive more realistic estimates of the amplitude uncertain- 
ties for the F03 data set. To this end, we adopted the recipe of Montgomery & O'Donoghue 
(1999), which suggests that the amplitude uncertainty should be equal to about 0.8 times 
the 1(7 noise level, irrespective of the actual amplitude of a mode. This leads to a{UV) — 
1.22 mmag, cr{B) — 1.22 mmag, cr(i?) = 0.90 mmag, and a{NIR) = 1.10 mmag. We note 
that these correspond to values 10 to 25 times larger than the formal fit errors quoted by 
F03 in their Tables 2 and 3, a huge difference perhaps in line with the reluctance shown by 
K98 to quote realistic errors for this type of data. While the approach of Montgomery & 
O'Donoghue (1999) is conservative, we believe that it leads to more realistic estimates of the 
amplitude uncertainties in this case, and we have consequently redone our analysis of the 
F03 data set with these revised values. 

Our results arc summarized in Table 5. Partial mode discrimination is possible for the 
three highest amplitude modes considered, and we show the corresponding fits in Figures 
5, 6, and 7. The best result is obtained for the largest amphtude mode observed by F03 
(481.75 s), for which, according to our analysis, the I value is equal to either 0, 1, or 2. 

3.4. HS 0039-1-4302 (J04) 

HS 0039+4302 has been observed by J04 using the 3-channcl ULTRACAM camera 
attached to the 4.2-m William Herschel Telescope. Some 16.2 h of u'g'r' Sloan photometry 
has been acquired on this relatively faint EC 14026 pulsator ai V — 15.5. Six modes with 
u'g'r' amplitudes have been isolated by the authors. While they have quoted an average 
amplitude uncertainty of 0.10 mmag for the 3 bandpasses of interest for all the 6 modes 
in their Table 3, we can use the results of JOS and Aerts et al. (2006) based on the same 
detector/telescope combination to infer wavelength-dependent amplitude uncertainties of 
(t{u') = 0.156 mmag, cr[g') = 0.063 mmag, and o'(r') = 0.081 mmag. These are the values 
that we used in our modeling of the 6 different modes. Note that during this modeling 
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exercise, we convolved the transmission curves of the Sloan filters it', g\ and r' with the 
response curves of the CCD's used in VhTKACAM kindly provided to us by V. Dhillon 
(2005, private communication). 

The results of our analysis for HS 0039+4302 are summarized in Table 6. They are 
interesting in that they show, despite the relative faintness of the target, that partial mode 
identification is possible. Except for the 134.44 s mode (the lowest amplitude one) for which 
mode discrimination is not feasible, partial identification is indeed possible. For instance. 
Figures 8, 9, and 10, corresponding to the three largest amplitude modes, indicate that these 
pulsations have values of the degree index / of either 0, 1, or 2. We note that the preliminary 
/ assignments made by J04 for the 6 modes they observed in HS 0039+4302 are consistent 
with our quantitative results. 



3.5. BAL 090100001 (O05) 

Oreiro et al. (2004) reported the discovery of short-period luminosity variations in 
BAL 090100001, the brightest and largest amplitude EC 14026 pulsator so far discovered. 
From a strictly observational point of view, these characteristics make it an ideal target 
to attempt multicolor photometry. Hence, less than a year after the initial discovery of the 
variabihty of BAL 090100001, O05 presented follow-up observations, including 9.3 h of simul- 
taneous Stromgren photometry on that star. This data set was gathered with the dedicated 
Stromgren photometer attached to the 0.9-m telescope at the Sierra Nevada Observatory. 

Our attempts to model the amplitude- wavelength behaviors of the two pulsation modes 
reported by O05 (see their Table 5) were frustrated and no acceptable fits were found for any 
of the assumed I values from to 5. Table 7 summarizes our results and clearly indicates 
that our quality-of-fit values Q are completely unacceptable. Figure 11 illustrates graphi- 
cally the situation for the largest amplitude mode (356.3 s). Again, as in the case of PC 
1605+072 discussed above, we cannot exclude completely the possibility that our modeling 
is inadequate for BAL 090100001, but we rather strongly suspect that the quoted amplitude 
uncertainties in Table 5 of O05 (most likely formal fit errors) underestimate largely the true 
errors. Unfortunately, unhke the case of PG 1605+072 in F03, not enough information is 
given in O05 that would have allowed us to obtain perhaps more realistic estimates of the 
amplitude uncertainties. Thus, we must conclude that we cannot model adequately the data 
of O05 as presented by them. 
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3.6. BAL 090100001 (BOS) 

BAL 090100001 was also observed by BOS who, in a major effort, gathered more than 
126 h of UBVR photometry on the 0.6-m telescope of the Mt. Suhora Observatory in Poland. 
The authors uncovered many periodicities, including not only short-period pulsations char- 
acteristic of EC 14026 stars, but also long-period oscillations most likely due to g-mode 
pulsations as found in the PG 1716 stars (and see also O05). We have concentrated our 
modehng effort on the 9 short-period modes (labelled /i through /g in BOS) with rehable 
amplitude determinations in all four bandpasses. We note that the authors provided a very 
careful frequency analysis and that they adopted the approach of Montgomery & O'Donoghue 
(1999) for estimating the uncertainties on the reported amplitudes and phases. For the am- 
plitudes, this leads to o"(f/) = 0.2S mmag, cr(i?) = 0.18 mmag, cr(V^) = 0.16 mmag, and (t{R) 
= 0.17 mmag, irrespective of the actual amplitudes of a given mode. 

The light curve of BAL 090100001, as observed by BOS, is dominated by a very large 
amplitude mode (the largest ever observed so far in a sdB pulsator) with a period of 3S6.19 
s. Its color amplitudes are U = 7S.23 mmag, B = S7.71 mmag, V = S3. 34 mmag, and R 
— SO. 26 mmag (see Table 1 of BOS). Not surprisingly, BOS also report the detection of the 
first and second harmonic of that mode, the first with an amplitude larger than 6 of the 9 
modes we retained for analysis. The dominant mode appears to be isolated, whereas the 
three next largest ones (354.20 s, 354.01 s, and 353.81 s) appear to form an almost perfectly 
symmetric triplet in frequency space, leading to the suggestion that this triplet could be 
due to rotational splitting. Considering that the main mode does not show an equivalent 
multiplet structure, this led BOS to suggest that the 3S6.19 s mode could be a / = mode, 
while the triplet would correspond to a rotationally-split / = 1 pulsation. It is obviously of 
high interest to verify if these sensible suggestions could be proven true. 

To model the data of BOS, we adopted the atmospheric parameters given in Table 2 (Tcs 
= 29,810 K and log g = S.S8), and we convolved our basic monochromatic quantities with 
the standard transmission curves of the Johnson/Cousins filters, the KPNO extinction curve, 
and a gray response for the CCD detector since, in the latter case, we had no information 
as to the exact response curve of the detector. The results of our effort for the 9 modes 
retained are summarized in Table 8. In addition, we include here Figures 12, 13, 14, and 15 
that refer, respectively, to the main mode and the triplet mentioned above. The other modes 
have relatively low amplitudes and it becomes increasingly difficult to discrimate between 
the possible values of / for them. 

Figure 12 illustrates our results for the main mode. Formally speaking, all model fits 
shown in the figure for that mode must be rejected since Q <^ 0.001 for all values of I 
considered in our simulations. At the same time, the figure also clearly illustrates that the 
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solution must be either / = or / = 1. We note that if the reported amphtude uncertainties 
are multiphed by a factor of 3.5 (to mimic the possibihty that these uncertaities have been 
perhaps underestimated), then the solutions / = and / = 1 become formally acceptable, 
while the others can be safely discarded. However, we have no particularly good reason to 
believe that the uncertainties on the amplitudes would have been underestimated by such a 
relatively large amount in this data set. Let us consider instead the data at face value and 
examine if we can rise up to the challenge offered by this very high S/N observation. 

The first obvious possibility to improve the match between the data points and one 
(/ = 0) or the other (/ = 1) of the possible solutions is to vary the model parameters within 
reasonable ranges, in particular the values of the atmospheric parameters T^g and log g that 
we assumed for BAL 090100001. The values of these parameters, as derived by Green et 
al. (2006), are T^s = 29,810 ± 400 K and log g = 5.58 ± 0.05. Exphcit calculations within 
these bounds indicate that the situation is not changed at the qualitative level: the values 
of for the I = and I — 1 solutions remain within a factor of 2 of each other as in Table 
8, and the values of Q are not significantly improved. 

The same is true by redoing the modeling exercise using, this time, the different set of 
atmospheric parameters derived by Oreiro et al. (2004), i.e., T^s — 29,450 ± 500 K and 
log g = 5.33 ± 0.10. In that case, the solution must again be either / = or / = 1, and 
the respective values are 1.39 x 10^ and 1.52 x 10^, both comparable to each other, but 
still way too high to be formally acceptable as their associated Q values arc much smaller 
than the criterion level of 0.001. We also experimented with changing the extinction curve 
from that of KPNO to one appropriate at sea level and found very little changes in our 
theoretical amplitudes. Likewise, assuming a specific CCD response, instead of using a fiat 
gray response as we did in our initial calculations, also led to very little qualitative changes 
compared to the situation depicted in Figure 12 and Table 8. (In that latter experiment, we 
did not use the response curve of the CCD employed by B05 since it was not available to 
us, but that of CCD21 of Steward Observatory as a surrogate.) Thus, we are unable to find 
a formally acceptable {Q > 0.001) / model for the main pulsation mode in BAL 090100001, 
but we still conclude that it must be either a / = or a / = 1 oscillation. 

This failure to model properly the main oscillation observed in BAL 090100001 is some- 
what bothersome (assuming again that the reported amplitude uncertainties are realistic 
estimates). We remind the reader that Randall et al. (2005), using the exact same tools 
as those used in this paper, have been able to model successfully the 182.42 s mode UL- 
TRACAM data reported by Jeffery et al. (2004) for KPD 2109+4401. We can think of a 
potentially important difference between the two cases, however, and it is the fact that the 
amphtude of the dominant mode in BAL 090100001 is so large that nonhnear effects are 
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obviously present in the form of the first and second harmonics. No such nonhnear features 
have been observed in KPD 2109+4401. The theory used to determine amphtude ratios 
between different bandpasses is strictly linear, and it is possible, although this remains un- 
proven, that it becomes inadequate to treat very large amplitude oscillations such as the 

356.19 s mode in BAL 090100001. 

Figures 13, 14. and 15 show our model fits for the triplet of modes which has been 
interpreted as a rotationally-split / = 1 mode by BOS. Taking again the data at face value, 
we find (sec Table 8) that the 354.20 s mode has / = 2, the 354.01 s mode has a value of / of 
cither 1 or 2, and the 353.81 s component has a value of / of cither 1, 2, or 4. If, however, 
the amplitude uncertainties are somewhat underestimated, then the / = 1 possibility for the 

354.20 s mode should not be dismissed too hastily because its Q value would not be much 
smaller than the passage value of Q — 0.001. Because of this, we prefer to be conservative 
and conclude that if the observed triplet of modes is indeed due to rotational splitting, then 
it must have a value of either ^ = 1 or / = 2. In the latter case, an unfavorable observation 
angle or some other cause could perhaps hide two components of the quintuplet. 

We have thus been unable to confirm the suggestion of B05 that the dominant 356.19 s 
mode in BAL 090100001 is a radial mode, and that the 354.20, 354.01, and 353.81 s modes 

are the components of a rotationally-split / = 1 mode. Our results are nevertheless consistent 
with that suggestion, but we found that the dominant mode could also have a degree index 
of / = 1 while the triplet could be three components of a split / = 2 mode. 

3.7. PG 0014+067 (JOS) 

The case of PG 0014+067 is particularly interesting because a specific asteroseismic 
model has been proposed by Brassard et al. (2001) including, of coTirsc, complete mode 
identification. This model has been refined subsequently by Charpinct ct al. (2005c) who 
exploited a higher quality (white light) data set and found the same mode identification. In 
both cases, in the search in parameter space for the optimal model, it was explicitly assumed 
that the observed modes had to belong to degree indices from / = up to and including 
I — 3. The inclusion of Z = 3 was deemed necessary because the observed mode density in 
PG 0014+067 is too large to be explained solely in terms of modes with / = 0, 1, and 2. In 
other words, there are more observed modes in the relevant period window of PG 0014+067 
than there are theoretical periods with / = 0, 1, and 2, so there seemed to be no choice but 
conclude that some additional modes with / > 3 are excited to visible levels in that star. 
Except for the restriction to / values from to 3, no a priori constraints were imposed in 
the asteroseismological exercises carried out by Brassard et al. (2001) and Charpinet et al. 
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(2005c), and the actual mode identification (the determination of both the radial order k 
and the degree index I for each observed pulsation) came out as a natural byproduct of the 
method. 

The results of Brassard et al. (2001) correspond to the first claimed successful astero- 
seismological exercise for a pulsating sdB star, and the importance of testing the proposed 
seismic model of PG 0014+067 cannot be overstated. It is in this spirit that JOS decided 
to attempt testing the / identification inferred in Brassard et al. (2001) using multicolor 
photometry. Even though the authors of JOS realized fully that, at ~ 16, PG 0014+067 
was going to be a challenge, they thought the matter to be sufficiently important to invest 
almost 30 h of WHT time using the ULTRACAM camera. It would have been difficult to do 
better than they did on this from an observational point of view. As in the previous targets 
studied in this paper, the details of their data set are given in Table 1. JOS isolated some 
10 distinct pulsation modes in PG 0014+067, and we modeled each one of these pulsations 
using the parameters hsted in Table 2. 

Not unexpectedly, our analysis reveals that most of the observed modes in PG 0014+067 
cannot be tested for I identification as the measured amplitudes usually do not have high 
enough S/N ratios. The details of our quantitative analysis are presented in Table 9 for the 
10 modes of interest. At the same time, our work also indicates that partial / discrimination 
is possible for the three largest amplitude modes uncovered by JOS, and we show our model 
fits for these modes in Figure 16 (141.01 s). Figure 17 (141.06 s), and Figure 18 (146. SO s). 
We find that the observed amplitudes of the 141.01 s mode can be accommodated if the 
degree index / has a value of either 0, 1, 2, or 4. This agrees with the conclusion of JOS that 
this mode (the /12 mode in their notation) cannot have a value of I — 3, and indeed we find 
that Q — 1.2x 10~^^ for that model, much too small to be acceptable. We also find that the 
141.06 s mode (/n in JOS) has a value of I of either 0, 1, or 2, again in agreement with the 
inference made in JOS. Finally, we find that the observed amplitude- wavelength distribution 
of the 146. SO s mode (/g in JOS) can be quantitatively explained if the mode has a I value 
of either 0, 1, or 2. This again is consistent with JOS who found that this mode cannot be a 
mode with / = 3 or / = 4. 

The most interesting result of this analysis is the conclusion, first put forward by JOS, 
that the 141.06 s mode (observed as the 141.07 s mode in the lower resolution white light 
data of Brassard et al. 2001) cannot have a degree index of / = 3. Given that Brassard et al. 
(2001) formally identified this mode as a / = 3 pulsation, this implies that a reevalution of 
their seismic model is warranted. It should be pointed out that at the time of the analysis 
of Brassard et al. (2001), the theory of multicolor photometry had not yet been apphed 
to realistic models of pulsating sdB stars. Using an Eddington limb darkening law (instead 
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of an exact form coming from detailed sdB model atmospheres), Brassard et al. (2001) 

estimated that the visibility factor would be equal to 1.0000, 0.7083, 0.3250, 0.0625, and 
0.0208 respectively for modes with / = 0, 1, 2, 3, and 4. Hence, the / = 3 modes would 
be more than three times "more visible" than their / = 4 counterparts, and it was deemed 
natural to limit the search to modes with / = 3 beyond those with / = 0, 1, and 2 since this 
was now sufficient to account for the relatively high density of observed modes. 

What we learned in Randall et al. (2005), among many other things, is that the visi- 
bility of a p-mode with / = 4 is, in fact, substantially larger than that of a / = 3 mode in 
a sdB pulsator in the optical domain. The discussion of Figure 1 in Randall et al. (2005) 
is quite explicit about this. Hence, it will be necessary in future asteroseismological exer- 
cises concerning EC 14026 pulsators (especially those showing more modes than could be 
accommodated by invoking the presence of only modes with / = 0, 1 , and 2) to include the 
possibility that the detected modes could also have values of Z = 4 as they have a higher 
probability of being detected than the / = 3 modes. In particular, a rcanalysis of the data 
presented by Brassard et al. (2001) and Charpinet et al. (2005c) on PG 0014+067 including 
this possibility must be carried out. This is being done and will be reported in due time. 

4. DISCUSSION 

We have presented in this paper a quantitative and homogeneous analysis of the broad- 
band multicolor photometric data sets gathered so far on rapidly pulsating sdB (EC 14026) 
stars. This analysis was carried out within the theoretical framework developed by Randall 
et al. (2005) with the goal of partial mode identification, i.e., the determination of the degree 
index / of each of the observed pulsation modes. With the exception of the very high S/N 
ratio u'g'r'/ULTRACAM data set obtained by Jeffery et al. (2004) on KPD 2109+4401 (an 
analysis of which will be presented in a separate publication), we considered all available 
data that we could find in the literature. This consists of 7 distinct data sets pertaining to 
6 different EC 14016 stars and involving 41 pulsation modes. 

Our final results are summarized in Table 10 where we list the name of the target star 
(1st column), the period of the mode of interest (2nd column), the acceptable values of / for 
that mode (3rd column), and a qualifier as to the level of mode discrimination achieved (4th 
column) . It can be seen that no completely convincing and unambiguous / identification has 
been possible on the basis of the available data, although we came close to that goal in the 
case of two modes in KPD 2109+4401 (K98) and three modes in BAL 090100001 (B05) for 
which only two possible values of I were found to be consistent with the observed amplitude- 
wavelength distributions. At the same time, we found that no mode discrimination (between 
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values from / = to Z = 5) was possible for 16 out of 41 modes. There is no doubt that 
these disappointing (but perhaps not totally unexpected) results are due to the fact that the 
available data sets were of insufficient sensitivity. The exception is the data set of BOS which 
features a remarkably high S/N value for the largest amplitude mode (356.19 s) uncovered 
in that study. We speculate here that nonlinear effects (observed in the form of the first 
and second harmonic of that mode) may have undermined our ability to model properly the 
amplitude-wavelength behavior of that particular mode. The latter shows an exceptionally 
large amphtude by EC 14026 star standards. We recall, in this context, that the feasibihty 
of our approach has been demonstrated beyond any doubt by Randall et al. (2005) in their 
analysis of the 182.42 s mode measured by Jeffery et al. (2004) in KPD 2109+4401. Hence, 
the case of the 356.19 s mode in BAL 090100001 remains enigmatic. 

It is interesting to examine the statistics of the numbers shown in Table 10. Considering 
only those cases where partial mode discrimination has been possible, we find that the I — 
solution comes up 19 times, the / = 1 solution 23 times, the / = 2 solution 24 times, the I = 3 
solution 3 times, the / = 4 solution 13 times, and the I — 5 solution once. While no single / 
index identification has been achieved in our study, one can argue, from a strictly statistical 
point of view, that these results are certainly consistent with the view that the pulsations 
seen in rapidly pulsating sdB stars arc low-degree modes, mostly with values of / = 0, 1, 
and 2. This is not a great surprise by any means, but it confirms the early interpretation of 
the EC 14026 phenomenon in terms of low-order, low-degree p-mode pulsations as presented 
by Charpinet et al. (1997). Moreover, the higher frequency of / = 4 solutions compared to 
that of Z = 3 solutions is consistent with the finding of Randall et al. (2005) that the former 
modes are more visible than the latter (in the optical domain). 

In conclusion, we find that in order to exploit the full potential of multicolor photom- 
etry for EC 14026 pulsators, it will be necessary in the future to gather higher S/N ratio 
observations than has been done up to now. This is certainly within reach, however, as our 
results clearly imply. Furthermore, it will be necessary in future asteroseismological exercises 
such as those carried out by Brassard et al. (2001) and Charpinet et al. (2005c) to include, 
as required, theoretical modes with / = 4 since those have a higher visibility factor in sdB 
stars than I = 3 modes. 

This work was supported in part by the NSERC of Canada. G. Fontaine also acknowl- 
edges the contribution of the Canada Research Chair Program. P. Bergeron is a Cottrell 
Scholar of the Research Corporation. 



Table 1. Basic characteristics of the data sets used in this paper 
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Table 2. Atmospheric parameters used in the computations of the theoretical amphtudes 
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Table 3. Fits of predicted UBVR amplitudes to those observed for the pulsation modes of 

KPD 2109+4401 reported by K98 
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Table 3 — Continued 
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0.709 


0.754 


0.834 


5 


1.85e+00 


6.03e- 


-01 


198.19 


4.800±0.105 


4.000±0.105 


4.100±0.085 


4.000±0.155 












5.318 


3.989 


3.733 


3.622 





4.90e+01 


1.33e- 


-10 




5.224 


4.021 


3.794 


3.694 


1 


3.32e+01 


2.91e- 


-07 




5.108 


4.055 


3.863 


3.765 


2 


1.90C+01 


2.77c- 


-04 




5.832 


3.771 


3.387 


2.902 


3 


2.22C+02 


7.82c- 


-48 




4.614 


4.177 


4.092 


4.087 


4 


6.29C+00 


9.85c- 


-02 




3.880 


4.073 


4.330 


4.795 


5 


l.lle+02 


7.10e- 


-24 



-23- 



-24- 



Table 4. Fits of predicted UBV amplitudes to those observed for the pulsation modes of 

PG 1618+563B reported by SOO 



Period 


U 


B 


V 


I 




Q 




(s) 


(mmag) 


(mmag) 


(mmag) 










139.30 


1.500±0.500 


1.200±0.600 


0.600±0.500 












1.292 


1.050 


0.983 





8.23e-01 


6.63e- 


-01 




1.276 


1.056 


0.995 


1 


8.82e-01 


6.43e- 


-01 




1.254 


1.062 


1.008 


2 


9.61e-01 


6.19e- 


-01 




1.442 


1.007 


0.853 


3 


3.74C-01 


8.30c- 


-01 




1.168 


1.079 


1.055 


4 


l.Slc+OO 


5.19c- 


-01 




1.078 


1.063 


1.105 


5 


1.79e+00 


4.09e- 


-01 


143.90 


4.600±0.700 


4.300±0.800 


2.600±1.000 












4.608 


3.745 


3.504 





1.30e+00 


5.23e- 


-01 




4.566 


3.778 


3.559 


1 


1.35e+00 


5.10e- 


-01 




4.508 


3.820 


3.627 


2 


1.43e+00 


4.89e- 


-01 




4.940 


3.478 


2.949 


3 


1.41e+00 


4.93c- 


-01 




4.278 


3.958 


3.872 


4 


2.01e+00 


3.66e- 


-01 




4.048 


3.994 


4.159 


5 


3.20e+00 


2.02e- 


-01 



-25- 



Table 5. Fits of predicted BUSCA amplitudes to those observed for the pulsation modes 

of PG 1605+072 reported by F03 



Period 


UV 


B 


R 


NIR 


I 




Q 




(s) 


(mmag) 


(mmag) 


(mmag) 


(mmag) 










440.51 


8.92dz 1.22 


6.85± 1.22 


6.85± 0.90 


6.66± 1.10 












8.80 


7.15 


6.99 


6.29 





2.10C-01 


9.76c- 


-01 




8.62 


7.18 


7.03 


6.41 


1 


2.28C-01 


9.73c- 


-01 




8.38 


7.22 


7.08 


6.51 


2 


3.72e-01 


9.46e- 


-01 




10.21 


7.55 


6.44 


3.11 


3 


1.21e+01 


6.97e- 


-03 




/ .00 


1 .10 


7 1 R 


7 10 
/ .iZ 


A 


i.Dze-|-uu 


D.ooe— 


-Ui 




6.04 


6.40 


6.99 


8.44 


5 


8.40e+00 


3.84e- 


-02 


475.82 


19.76d= 1.22 


16.94± 1.22 


16.68± 0.90 


17.18± 1.10 












21.21 


17.19 


16.80 


15.10 





5.09C+00 


1.66c- 


-01 




20.73 


17.27 


16.92 


15.42 


1 


3.36e+00 


3.40e- 


-01 




20.17 


17.41 


17.07 


15.70 


2 


2.26e+00 


5.19e- 


-01 




24.16 


18.05 


15.40 


7.58 


3 


9.26e+01 


5.92e- 


-20 




1 8 on 
io.zu 


1 7 08 
1 / .Zo 


1 7 "^0 


1 7 07 
1 / .Z / 


A 


z.zoe-ruu 


o.zue— 


-Ui 




14.58 


15.52 


17.00 


20.62 


5 


2.95e+01 


1.79e- 


-06 


481.75 


36.88± 1.22 


29.27± 1.22 


28.45± 0.90 


28.72± 1.10 












37.05 


30.03 


29.33 


26.36 





6.00e+00 


1.12e- 


-01 




36.20 


30.14 


29.53 


26.92 


1 


4.98C+00 


1.74c- 


-01 




35.16 


30.36 


29.77 


27.40 


2 


6.44C+00 


9.22c- 


-02 




42.50 


31.80 


27.13 


13.40 


3 


2.23e+02 


4.30e- 


-48 




31.62 


30.06 


30.13 


30.06 


4 


2.41e+01 


2.34e- 


-05 




25.18 


26.83 


29.41 


35.70 


5 


1.38e+02 


9.56e- 


-30 


503.70 


22.59± 1.22 


22.45± 1.22 


18.79± 0.90 


18.02± 1.10 












24.48 


19.80 


19.34 


17.36 





7.83e+00 


4.96c- 


-02 




23.89 


19.89 


19.48 


17.76 


1 


6.18C+00 


1.03c- 


-01 




23.18 


20.05 


19.67 


18.11 


2 


5.04e+00 


1.69c- 


-01 




28.05 


21.12 


18.03 


9.00 


3 


8.97e+01 


2.53c- 


-19 




20.84 


19.86 


19.92 


19.90 


4 


l.lle+01 


1.12e- 


-02 



-26- 



Table 5 — Continued 



Period 


UV 


B 


R 


NIR 


I 




Q 






(mmag) 


(mmag) 


(mma.g) 


(mmag) 












16.52 


17.65 


19.39 


23.61 


5 


6.68e+01 


2.09e- 


-14 


528.70 


7.65± 1.22 


8.25± 1.22 


6.79± 0.90 


6.88± 1.10 












8.84 


7.13 


6.96 


6.24 





2.16e+00 


5.39e- 


-01 




8.62 


7.17 


7.03 


6.41 


1 


1.66e+00 


6.45e- 


-01 




8.35 


7.24 


7.10 


6.55 


2 


1.22C+00 


7.48c- 


01 




9.99 


7.58 


6.47 


3.27 


3 


1.50C+01 


1.83c- 


03 




7.51 


7.19 


7.21 


7.22 


4 


1.09C+00 


7.80c- 


-01 




5.96 


6.40 


7.05 


8.61 


5 


6.79e+00 


7.91e- 


-02 



-27- 



Table 6. Fits of predicted ULTRACAM amplitudes to those observed for the pulsation 

modes of HS 0039+4302 reported by J04 



Period 


u' 


g' 


r' 


I 




Q 




(s) 


(mmag) 


(mmag) 


(mmag) 










134.44 


0.890±0.156 


0.730±0.063 


0.600±0.081 












0.914 


0.690 


0.660 





9.81C-01 


6.12c- 


-01 




0.900 


0.690 


0.663 


1 


9.99C-01 


6.07c- 


-01 




0.886 


0.693 


0.666 


2 


1.02e+00 


6.02e- 


-01 




1.178 


0.676 


0.545 


3 


4.62e+00 


9.95e- 


-02 




n 7Q/I 
u. / y4 


u.oy / 


U.Do / 


A 


1 sio_i_nn 
i.oie-|-uu 


^.U'ie— 


-Ui 




0.678 


0.676 


0.735 


5 


5.38e+00 


6.80e- 


-02 


180.13 


1.590±0.156 


1.350±0.063 


1.410±0.081 












1.820 


1.363 


1.300 





4.06C+00 


1.31c- 


-01 




1.780 


1.367 


1.311 


1 


3.06e+00 


2.16e- 


-01 




1.732 


1.374 


1.321 


2 


2.18e+00 


3.36e- 


-01 




2.152 


1.358 


1.093 


3 


2.83e+01 


7.12e- 


-07 






l.ooo 


1 '^70 
l.o / Z 


A 


D.u^te — ui 


7 '^Qo 

/ .oye— 


-ui 




1.322 


1.340 


1.487 


5 


3.89e+00 


1.43e- 


-01 


181.89 


3.590±0.156 


2.960±0.063 


2.640±0.081 












3.820 


2.859 


2.728 





5.90e+00 


5.23e- 


-02 




3.734 


2.868 


2.749 


1 


4.82e+00 


8.98c- 


-02 




3.628 


2.878 


2.768 


2 


4.22e+00 


1.21c- 


-01 




4.542 


2.874 


2.313 


3 


5.54eH-01 


9.15e- 


-13 




3.216 


2.890 


2.867 


4 


1.48e+01 


6.10e- 


-04 




2.750 


2.788 


3.097 


5 


6.83e+01 


1.51e- 


-15 


182.79 


2.390±0.156 


2.000±0.063 


1.800±0.081 












2.582 


1.932 


1.843 





2.95C+00 


2.28c- 


-01 




2.524 


1.938 


1.858 


1 


2.21e+00 


3.31c- 


-01 




2.452 


1.946 


1.871 


2 


1.67e+00 


4.35c- 


-01 




3.064 


1.941 


1.562 


3 


2.82e+01 


7.66c- 


-07 




2.174 


1.954 


1.939 


4 


5.37e+00 


6.82e- 


-02 



-28- 



Table 6 — Continued 



Period 


u' 


9' 


r' 


I 




Q 




(s) 


(mmag) 


(mmag) 


(mmag) 












1.860 


1.886 


2.096 


5 


2.81e+01 


7.78e- 


-07 


192.58 


5.730±0.156 


4.350±0.063 


4.050±0.081 












5.764 


4.309 


4.108 





9.90C-01 


6.10c- 


-01 




5.624 


4.320 


4.140 


1 


1.93e+00 


3.80e- 


-01 




5.454 


4.337 


4.170 


2 


5.38e+00 


6.79e- 


-02 




6.776 


4.351 


3.501 


3 


9.096+01 


1.79e- 


-20 




4.824 


4.352 


4.321 


4 


4.49e+01 


1.77e- 


-10 




4.122 


4.191 


4.671 


5 


1.71e+02 


6.19e- 


-38 


234.11 


6.110±0.156 


4.930±0.063 


4.620±0.081 












6.536 


4.843 


4.607 





9.38e+00 


9.17e- 


-03 




6.328 


4.862 


4.658 


1 


3.32C+00 


1.90e- 


-01 




6.072 


4.888 


4.704 


2 


1.59e+00 


4.53e- 


-01 




7.290 


4.953 


3.993 


3 


1.17e+02 


3.37e- 


-26 




5.336 


4.904 


4.890 


4 


3.59e+01 


1.60e- 


-08 




4.554 


4.699 


5.315 


5 


1.87e+02 


3.07e- 


-41 



-29- 



Table 7. Fits of predicted uvby amplitudes to those observed for the pulsation modes of 

BAL 090100001 reported by O05 



Period 


u 


V 


b 


y 






Q 




(mmag) 


(mmag) 


(mmag) 


(mmag) 








351.70 


16.73± 0.08 


17.43± 0.09 


16.84± 0.09 


15.76± 0.08 










21.02 


13.97 


15.27 


13.27 





5.62e+03 


O.OOe+OO 




20.24 


14.63 


15.63 


14.07 


1 


3.51e+03 


O.OOe+00 




19.37 


15.30 


15.93 


14.83 


2 


1.89e+03 


O.OOe+00 




20.87 


14.83 


15.03 


13.15 


3 


4.99C+03 


O.OOc+00 




17.27 


16.26 


16.60 


16.32 


4 


2.71C+02 


1.90C-58 




14.49 


16.21 


16.29 


18.70 


5 


2.35e+03 


O.OOe+00 


356.30 


49.36± 0.11 


41.34± 0.11 


38.55± 0.10 


38.03± 0.10 










54.58 


36.23 


39.63 


34.40 





5.84e+03 


O.OOe+00 




52.05 


37.63 


40.20 


36.19 


1 


2.34e+03 


O.OOe+00 




49.34 


39.00 


40.62 


37.81 


2 


8.86e+02 


8.47-192 




54.02 


38.46 


38.99 


34.10 


3 


4.04e+03 


O.OOe+00 




43.19 


40.71 


41.56 


40.88 


4 


4.90e+03 


O.OOe+00 




35.77 


40.04 


40.20 


46.21 


5 


2.24e+04 


O.OOe+00 



-30- 



Table 8. Fits of predicted UBVR amplitudes to those observed for the pulsation modes of 

BAL 090100001 reported by BOS 



Period 


U 


B 


V 


R 






Q 


(s) 


(mmag) 


(mmag) 


(mmag) 


(mmag) 








356.19 


75.23± 0.25 


57.71± 0.18 


53.34± 0.16 


50.26± 0.17 










77.92 


58.14 


52.05 


49.37 





2.15e+02 


2.68C-46 




73.18 


57.96 


53.38 


51.31 


1 


1.07C+02 


3.84C-23 




68.50 


57.76 


54.53 


52.85 


2 


l.Ole+03 


2.91-219 




77.85 


59.92 


52.99 


46.11 


3 


8.61e+02 


2.11-186 




Oo.oo 


OD.OO 


00. ys 


i^R 07 


A 


0.0/ e-|-Uo 


u.uue-|-uu 




46.98 


49.55 


55.85 


63.17 


5 


2.08e+04 


O.OOe+00 


354.20 


26.71± 0.25 


21.92± 0.18 


20.53± 0.16 


19.92± 0.17 










29.52 


22.03 


19.73 


18.71 





2.02C+02 


1.57e-43 




27.77 


22.00 


20.26 


19.47 


1 


2.81C+01 


3.55e-06 




26.04 


21.94 


20.72 


20.08 


2 


9.53e+00 


2.30e-02 




29.50 


22.69 


20.06 


17.46 


3 


3.61e+02 


7.16e-78 






01 f^f; 
Zl.OO 


01 'VX 

Zl.OO 


01 4/L 


A 


A (T^a-LnO 

4.uoe-rUZ 


D.lOe — o / 




17.98 


18.96 


21.37 


24.16 


5 


2.14e+03 


O.OOe+00 


354.01 


15.82± 0.25 


12.86± 0.18 


11.62± 0.16 


11.57± 0.17 










17.14 


12.79 


11.46 


10.87 





4.61C+01 


5.27e-10 




16.12 


12.77 


11.76 


11.30 


1 


4.91C+00 


1.78C-01 




15.11 


12.73 


12.02 


11.65 


2 


1.51C+01 


1.71e-03 




17.13 


13.17 


11.65 


10.14 


3 


l.Ole+02 


7.47e-22 




13.00 


12.49 


12.36 


12.43 


4 


1.78e+02 


1.91e-38 




10.41 


10.98 


12.38 


13.99 


5 


8.02e+02 


1.76-173 


353.81 


5.86± 0.25 


5.12± 0.18 


4.71± 0.16 


4.59± 0.17 










6.75 


5.04 


4.51 


4.28 





1.77e+01 


4.99C-04 




6.35 


5.03 


4.63 


4.45 


1 


4.99C+00 


1.73C-01 




5.96 


5.02 


4.74 


4.59 


2 


4.92C-01 


9.21C-01 




6.74 


5.18 


4.58 


3.99 


3 


2.57e+01 


1.13e-05 




5.13 


4.93 


4.88 


4.91 


4 


1.42e+01 


2.69e-03 
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Table 8 — Continued 



Period 


U 


B 


V 


R 




9 

X 


Q 




(s) 


(mmagj 


(mmagj 


(mmagj 


(mmagj 












4 1 9 


4 35 


4 90 


5 54 




q qQp_i_ni 


9 93p- 


-21 

Zf X 


35(1 3Q 




1 81+ n 18 


2 08+ 16 


1 45+ f) 1 7 












2 53 


1 89 


1 69 


1 61 

X . \J X 


n 


8 33p+00 


3 96r- 


-02 




2 38 


1 89 


1 74 

X . 1 t: 


1 67 

X . U 1 


1 

X 


\j . X WW 




-09 




9 9'^ 


1 88 


1 78 
X . 1 


1 79 


9 


u .0 xt.n^ww 


9 75p- 


-09 
wz 




2.54 


1.95 


1.72 


1.50 


3 


7.17e+00 


6.66e- 


-02 




1.93 


1.85 


1.83 


1.84 


4 


9.20e+00 


2.67e- 


-02 




1 54 


1 69 


1 83 

X .00 


9 07 


e; 

U 




2 14p- 

Zf . XrtC 


-05 




2 nn- n 25 


1 72+ n 18 

X.I Zj I L/ . XL? 


1 oil A 1C 
X.OJ 1— L/.XW 


1 40+ f) 17 

X .'-±0_l_ L/.XI 












2 15 


1 60 

X . \J\J 


1 44 

X . J. J. 


1 36 


n 


1 a2p+no 


6 1 1p- 

W. X xc 


-01 

yj X 




2 02 


1 60 

X . \J\J 


1 47 

X .'-± 1 


1 41 

X .rt X 


1 

X 




6 53p- 


-01 

yj X 






1 5q 


1 50 


1 4(S 


9 


9 9ipj_no 


5 '?1p- 


-01 

U -L 




9 1 5 




X . TlVj 


1 97 






4 1 5p- 


-01 

W X 




1.62 


1.56 


1.54 


1.55 


4 


5.48e+00 


1.40P- 


-01 




1.30 


1.37 


1.54 


1.74 


5 


1.64e+01 


9.45e- 


-04 


34Q 83 


1 04+ f) 25 


1 45+ f) ic 


1 50+ f) -|« 


1 45+ f) 1 7 












2 1 3 


1 59 


1 43 

X .'-±(J 


1 35 


n 

yj 


2 45p+no 


4 84p- 


-01 

yj X 




2 nn 


1 59 


1 46 


1 41 


1 

X 


1 26p+00 


7 38p- 


-01 




1 QQ 
i.OO 


i.Oo 


1 /in 
i.4y 


1 /I c; 
i.40 



Z 


o.yie — Ui 


o.zoc— 


ni 
-Ui 




2.13 


1.64 


1.45 


1.26 


3 


3.61e+00 


3.07P- 


-01 




1.62 


1.56 


1.54 


1.55 


4 


2.37C+00 


5.00c- 


-01 




1.30 


1.38 


1.55 


1.75 


5 


9.80e+00 


2.04e- 


-02 


264.82 


5.94± 0.25 


4.73± 0.18 


4.16± 0.16 


4.07± 0.17 












6.09 


4.64 


4.21 


4.01 





8.73e-01 


8.32e- 


-01 




5.84 


4.62 


4.27 


4.10 


1 


l.OOe+00 


8.01e- 


-01 




5.59 


4.62 


4.34 


4.19 


2 


4.15e+00 


2.46e- 


-01 




6.45 


4.73 


4.17 


3.64 


3 


1.04e+01 


1.55e- 


-02 
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Table 8 — Continued 



Period 


U 


B 


V 


R 






Q 




(s) 


(mmag) 


(mmag) 


(mmag) 


(mmag) 












4.83 


4.54 


4.46 


4.46 


4 


2.96e+01 


1.70e- 


-06 




3.85 


4.05 


4.47 


4.99 


5 


1.17e+02 


4.05e- 


-25 


264.08 


2.58± 0.25 


1.83± 0.18 


1.51± 0.16 


1.60± 0.17 












2.40 


1.82 


1.65 


1.57 





1.37e+00 


7.12e- 


-01 




2.29 


1.82 


1.68 


1.61 


1 


2.40e+00 


4.95e- 


-01 




2.19 


1.81 


1.70 


1.64 


2 


3.92e+00 


2.71e- 


-01 




2.54 


1.86 


1.64 


1.43 


3 


1.69e+00 


6.39e- 


-01 




1.89 


1.77 


1.74 


1.74 


4 


1.06e+01 


1.43e- 


-02 




1.50 


1.58 


1.74 


1.94 


5 


2.67e+01 


6.83e- 


-06 
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Table 9. Fits of predicted ULTRACAM amplitudes to those observed for the pulsation 

modes of PG 0014+067 reported by JOS 



Period 


u' 


g' 


r' 


I 




Q 




(s) 

V / 


(mmag) 


(mmag) 


(mmag) 










100.29 


0.900±0.200 


0.690±0.070 


0.600±0.100 












0.864 


0.674 


0.648 





3.20C-01 


8.52c- 


-01 




0.854 


0.674 


0.650 


1 


3.55C-01 


8.37c- 


-01 




0.842 


0.675 


0.650 


2 


3.87e-01 


8.24e- 


-01 




1.200 


0.662 


0.497 


3 


3.47e+00 


1.77e- 


-01 




U. / DD 


U.D/ D 


U.DD/ 


A 

4 


y.ooe — ui 


D.z / e— 


ni 

-Ui 




0.676 


0.661 


0.710 


5 


2.63e+00 


2.68e- 


-01 


139.14 


0.800±0.200 


0.710±0.080 


0.800±0.100 












0.952 


0.736 


0.707 





1.56C+00 


4.59c- 


-01 




0.934 


0.738 


0.711 


1 


1.36C+00 


5.05c- 


-01 




0.912 


0.741 


0.715 


2 


1.19e+00 


5.51e- 


-01 




1.162 


0.734 


0.565 


3 


8.89e+00 


1.18e- 


-02 




u.ozz 


n 7AA 


U. / oo 


A 
4 


0. / / e — ui 


7 /IQo 


-ui 




0.728 


0.725 


0.794 


5 


1.69e-01 


9.19e- 


-01 


140.98 


1.200±0.300 


l.lOOiO.lOO 


l.OOOiO.lOO 












1.370 


1.060 


1.017 





5.13e-01 


7.74e- 


-01 




1.342 


1.061 


1.022 


1 


4.25e-01 


8.08c- 


-01 




1.306 


1.062 


1.024 


2 


3.29e-01 


8.48c- 


-01 




1.754 


1.113 


0.857 


3 


5.47eH-00 


6.50e- 


-02 




1.166 


1.056 


1.048 


4 


4.34e-01 


8.05e- 


-01 




1.010 


1.007 


1.103 


5 


2.33e+00 


3.12e- 


-01 


141.01 


5.900±0.300 


5.030±0.090 


4.600±0.100 












6.336 


4.900 


4.702 





5.24C+00 


7.28c- 


-02 




6.202 


4.902 


4.721 


1 


4.51C+00 


1.05c- 


-01 




6.038 


4.909 


4.735 


2 


3.83C+00 


1.47c- 


-01 




8.078 


5.127 


3.948 


3 


9.63C+01 


1.20c- 


-21 




5.390 


4.883 


4.845 


4 


1.16e+01 


3.10e- 


-03 
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Table 9 — Continued 



Period 

(s) 


(mmagj 


9 

(mmagj 


r 

(mmagj 


I 


9 

X 


Q 






4.690 


4.675 


5.122 


5 


5.91e+01 


1.50e- 


-13 


141 nfi 


4 c;no+o 9no 


Q onn-Lf) nyn 
o.^yjyj^Liyj .yj 1 u 


o 040+0 OQO 














3.222 


3 092 





3 22r+00 


2 OOp- 


-01 

W -L 




4 089 


996 


S 1 08 


1 

± 


c; 07p-I-00 


7 Ql p- 


-09 




S Q7fi 


3 933 


3 118 


9 


7 83p+00 


9 OOp- 

Zi . w W^ 


-09 

W ^ 






S 9Q8 


9 540 






1 74p- 


-10 




3.564 


3.229 


3.204 


4 


2.54e+01 


3.07e- 


-06 




3.128 


3.118 


3.416 


5 


6.59e+01 


4.90e- 


-15 


1 4fi 50 


4 nno+0 200 


Kj .\Jkj\J -i—\J . \J\J\J 


2 fiSO+0 080 












3 834 


1 963 


2 842 





fi 05p+00 


4 85p- 


-02 










1 

± 






-09 




S 646 


9 Q6Q 


9 86'^ 


At 


q 4Qp_L00 


8 Q7p- 






4 766 


'\ 060 


9 S61 


Q 


Q 08p4-01 


9 04p- 


-07 




3.258 


2.959 


2.938 


4 


2.55C+01 


2.86e- 


-06 




2.858 


2.854 


3.133 


5 


7.32e+01 


1.25e- 


-16 


1 "sn 47 




f) fion+n 07n 


Koo+O 090 












876 


677 


649 





1 c;4p_01 


9 96p- 


-01 

yj ±. 




0.858 


0.678 


0.653 


-1 

1 


2.08e— Oi 


9.0ic- 


-01 




0.834 


0.680 


0.656 


2 


2.87C-01 


8.66e- 


-01 




1.076 


0.696 


0.538 


3 


1.62e+00 


4.45c- 


-01 




0.746 


0.679 


0.674 


4 


9.79e-01 


6.13c- 


-01 




0.654 


0.654 


0.719 


5 


2.98eH-00 


2.25e- 


-01 


150.78 


0.800±0.100 


0.570±0.070 


0.600±0.100 












0.774 


0.598 


0.573 





2.97e-01 


8.62e- 


-01 




0.762 


0.602 


0.580 


1 


3.99e-01 


8.19e- 


-01 




0.746 


0.608 


0.587 


2 


6.06e-01 


7.39e- 


-01 
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Table 9 — Continued 



reriod 


u 


9 




7 


2 


Q 


(s) 


(mmagj 


(mmagj 


(mmagj 










880 


570 


440 


3 


3 20e+00 

KJ • £^\J\J\^ \J\J 


2 n2p— 01 




684 


0.622 


618 


4 




■ J .owe \J X 




61 8 


61 8 


679 




A 4ip_i_nn 


1 IDp— 01 


154.94 


0.700±0.200 


0.670±0.060 


0.780±0.080 










0.910 


0.703 


0.674 





3.16e+00 


2.05e-01 




0.890 


0.704 


0.678 


1 


2.86e+00 


2.39e-01 




0.864 


0.705 


0.680 


2 


2.58e+00 


2.76e-01 




1.096 


0.714 


0.552 


3 


1.26e+01 


1.87e-03 




0.776 


0.707 


0.702 


4 


1.47e+00 


4.81e-01 




0.686 


0.687 


0.756 


5 


1.75e-01 


9.16e-01 


168.80 


0.600±0.200 


0.410±0.070 


0.400±0.100 










0.542 


0.418 


0.401 





9.74C-02 


9.52e-01 




0.530 


0.419 


0.403 


1 


1.41e-01 


9.32e-01 




0.514 


0.420 


0.405 


2 


2.10e-01 


9.00e-01 




0.644 


0.427 


0.331 


3 


5.85e-01 


7.46e-01 




0.460 


0.421 


0.418 


4 


5.47e-01 


7.61e-01 




0. 106 


0.108 


0.150 


5 


1.19e+00 


5.52e-01 
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Table 10. Summary of partial mode identification in pulsating sdB stars 



Object 


Period 


Possible / identification 


Mode discrimination 






(0-5) 




KPD 2109+4401 (K98) 


182.42 


0,1,2 


partial 




184.72 


0,1,2,4 


partial 






9 4 


T~\ r\ x'-t- 1 r\ 1 

yaL 11 di 




1Q1 R'^ 


1 9 4 


T~\ O T*-l- 1 O 1 




IcfKJ.O L 


n 1 9 


^dl tldi 




±cf\J.\Jcf 


n 1 9 4 ^ 






198.19 


2,4 


partial 


rKji ioio-roooij (^ouuj 


ioy.ou 




no 






n 1 9 'i 4 


IIU 






n 1 9 4 


IIU 






n 1 9 4 


pdl 11 di 




481.75 


0,1,2 


partial 




503.70 


0,1,2,4 


partial 




^98 7n 


n 1 9 4 ^ 


LiU 


HS 0039+4302 (J04) 


134.44 


0,1,2,3,4,5 


no 




180.13 


0,1,2,4,5 


partial 




181.89 


0,1,2 


partial 




182.79 


0,1,2,4 


partial 




192.58 


0,1,2 


partial 




234.11 


0,1,2 


partial 


BAL 090100001 (O05) 


351.70 


no acceptable fit 


no 




356.30 


no acceptable fit 


no 



BAL 090100001 (B05) 264.08 

264.82 
349.83 



0,1,2,3,4 
0,1,2,3 
0,1,2,3,4,5 



partial 
partial 
no 



Table 10 — Continued 



Object Period Possible / identification Mode discrimination 

(s) (0-5) 



350.18 


0,1,2,3,4 


partial 


350.39 


0,1,2,3,4 


partial 


353.81 


1,2,4 


partial 


354.01 


1,2 


partial 


354.20 


1,2 


partial 


356.19 


0,1 


partial 


PG 0014+067 (J05) 100.29 


0,1,2,3,4,5 


no 


139.14 


0,1,2,3,4,5 


no 


140.98 


0,1,2,3,4,5 


no 


141.01 


0,1,2,4 


partial 


141.06 


0,1,2 


partial 


146.50 


0,1,2 


partial 


150.47 


0,1,2,3,4,5 


no 


150.78 


0,1,2,3,4,5 


no 


154.94 


0,1,2,3,4,5 


no 


168.80 


0,1,2,3,4,5 


no 
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FIGURE CAPTIONS 

Fig. 1 — Fits to the U, B, V, and R pulsational amplitudes observed for the 182.42 s mode 
of KPD 2109+4401 by K98. The predicted amphtude-wavelength behaviors of modes with 
I — to I — 5 have been fit to the observed values using a least-squares procedure. The 
curves with I — 0, I — 1, and I — 2 provide acceptable fits according to the quality-of-fit Q 
quantity (see text). In contrast, the curves with I — 3, I — A, and / = 5 do not provide viable 
fits and, consequently, these possible I identifications must be rejected for that mode. 

Fig. 2 — Similar to Fig. 1, but for the mode with a period of 196.31 s. Here the acceptable 
fits are again for values of the degree index / = 0, Z = 1, and 1 — 2. 

Fig. 3 — Similar to Fig. 1, but for the mode with a period of 198.19 s. Here the only 
acceptable fit from a formal point of view is the one corresponding to / = 4. 

Fig. 4 — Fits to the U, B, and V pulsational amplitudes observed for the 139.3 s mode of 
PG 1618+563B by SOO. No mode discrimination is possible here. 

Fig. 5 — Fits to the BUSCA {UV, B, R, and NIR) pulsational amplitudes observed for the 
475.82 s mode of PG 1605+072 by F03. Partial mode discrimination is possible here, with 
the values / = 3 and / = 5 excluded. 

Fig. 6 — Similar to Fig. 5, but for the mode with a period of 481.75 s. Here the acceptable 
fits are for values of the degree index / = 0, / = 1, and 1 = 2. 

Fig. 7 — Similar to Fig. 5, but for the mode with a period of 503.70 s. Partial mode 
discrimination is possible, with the values / = 3 and I = 5 excluded. 

Fig. 8 — Fits to the ULTRACAM [u' , g' , and r') pulsational amphtudes observed for the 
181.89 s mode of HS 0039+4302 by J04. The acceptable values of I for that mode are I = 0, 
Z = 1, and / = 2. 

Fig. 9 — Similar to Fig. 8, but for the mode with a period of 192.58 s. Here the acceptable 
fits are for values of the degree index / = 0, Z = 1, and 1 = 2. 

Fig. 10 — Similar to Fig. 8, but for the mode with a period of 243.11 s. The acceptable fits 
are for values of the degree index / = 0, / = 1, and 1 = 2. 

Fig. 11 — Fits to the -u, h and y pulsational amplitudes observed for the 356.3 s mode of 
BAL 090100001 by O05. No mode discrimination is possible here. 

Fig. 12 — Fits to the [/, S, V ^ and R pulsational amplitudes observed for the 356.19 s mode 
of BAL 090100001 by B05. There is no acceptable formal fit, but it is clear that the mode 
must have a degree index / = or / = 1. 
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Fig. 13 — Similar to Fig. 12, but for the mode with a period of 354.20 s. Formally, the only 
acceptable {Q > 0.001) fit is that with / = 2, although the solution with I — 1 should not be 
discarded since the data points fall in between those two model curves. 

Fig. 14 — Similar to Fig. 12, but for the mode with a period of 354.01 s. The acceptable 
fits are for values of the degree index I — 1, and 1 — 2. 

Fig. 15 — Similar to Fig. 12, but for the mode with a period of 353.81 s. The acceptable 
fits are for values of the degree index I — 1, 1 — 2, and I — A. 

Fig. 16 — Fits to the ULTRACAM {u', g', and r') pulsational amplitudes observed for the 
141.01 s mode of PG 0014+067 by J05. The acceptable values of / for that mode are / = 0, 
i = 1, / = 2, and / = 4. 

Fig. 17 — Similar to Fig. 16, but for the mode with a period of 141.06 s. Here the acceptable 
fits are for values of the degree index / = 0, / = 1, and 1 = 2. 

Fig. 18 — Similar to Fig. 16, but for the mode with a period of 146.50 s. The acceptable 
fits are for values of the degree index / = 0, / = 1, and 1 = 2. 
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Figure 1 
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Figure 2 
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Figure 3 
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Figure 8 
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Figure 9 
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Figure 10 
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Figure 13 
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Figure 15 
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Figure 16 



-58- 




Figure 17 
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Figure 18 



